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FOREWARD 


A gas-surface Energy Transfer Experiment was flown on OGO-6 to 
measure the kinetic-energy flux of the upper atmosphere relative to the 
orbiting satellite. 2 Energy transfer between 10 microwatts/ cm^ and 
0.1 W/cm2 was measured by short-term frequency changes of temperature- 
sensitive quartz crystals used in the energy transfer probe. ^ The condition 
of the surfaces was continuously monitored by a quartz crystal microbalance 
to determine the effect surface contamination had on energy accommodation. 

Results from the experiment provided the first direct measurements 
of upper atmospheric interactions with surf aces, 4 as well as, new data on 
surface contamination. 5 Measurements showed that the energy transfer by 
the iqjper atmosphere to surfaces could be readily measured below 500 km. 

The energy flux varied greatly from day to night and with altitude and 
reached a maximum value of 110 microwatts/cm^ at the 400 km perigee of 
OGO-6. Contamination outgassing primarily from the solar panels signifi- 
cantly modified the aluminum surfaces. ° Theoretical analysis had pre- 
dicted energy- accommodation coefficients for aluminum of 0.9 but contami- 
nation limited this value to 0.2. 

The knowledge gained from the experiment has lead to two new instruments 
for space measurements. One instrument is the Upper Atmospheric Density 
Gauge. ^ The other instrument is the Theimoelectrically-Cooled Quartz Crystal 
Microbalance for measuring surface contamination as a function of tenqjerature.® 

Much of the data obtained from the experiment has been published.^"' 

This Final Report covers analysis not previously reported and new work 
completed this year. The Report is divided into three parts. 

Part 1 is entitled "Gas -Surface Interaction Studies". Results are 
given on the computer analysis and laboratory tests performed to optimize 
the operation of the energy transfer probe flown on OGO-6. 

Part 2 is entitled "Bombardment of OGO-6 Surfaces by High-Energy Par- 
ticles". Data are given on the penetration and trapping of protons in 
surfaces at energies up to an MeV. 

Part 3 is a paper entitled "Thermoelectrically-Cooled Quartz Crystal 
Microbalance" recently presented at the 7th Conference on Space Simulation, 

Los Angeles, California, November 1973. This new microbalance is to be 
used to monitor contamination on the Space Shuttle. Its development was 
partially funded by the NASA Marshall Space Flight Center under Contract 
NAS8-27879. 
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PART I 


GAS-SURFACE INTERACTION STUDIES 


D. McKeown, J.M, Bowyer, Jr, and R.S, Dummer 



GAS-SURFACE INTERACTION STUDIES 
D. McKeown, J. M. Bowyer, and R. S. Dummer 


1.0 INTRODUCTION 


When a molecular beam is reflected from a surface, it gives up part 
of its kinetic energy as a result of the interaction. At energies below a 
few hundred eV, most of the energy transferred is released as heat since atoms 
ejected from the surface by the beam will carry little of the energy^»2 away. 

In studying molecular beam interactions it is important to know what frac- 
tion of the energy is transferred. If the amount of energy transferred can 
be measured, the accommodation coefficient of the surface can be determined 
and the nature of the interaction better understood. 

The usual way of measuring energy transferred is by attaching a ther- 
mocouple to the target. This approach works well at high beam energy but at 
energies below one keV it does not have the required sensitivity. We have 
found as a result of working with quartz crystals to measure sputtering^ that 
energy transfer can also be detected by a crystal because its oscillating fre- 
quency is temperature dependent.'* By using special cuts of quartz, tempera- 
ture changes of 10 can be measured. 

There are several advantages in using quartz crystals to measure energy 
transfer. Quartz dissipates heat poorly and a small heat flux will produce a 
relatively large increase in the temperature of a crystal. Many different 
materials can be plated or bonded onto quartz. It is a stable substance and 
can be used in a vacuum at high tenqieratures . Its frequency -temperature coeffi- 
cient remains relatively constant down to O^K. Experiments are performed by 
placing a crystal in a molecular beam and measuring the energy transferred to 
it by monitoring changes in its oscillating frequency.^ 


2.0 


ENERGY TRANSFER 


The rate of change of the crystal tenperature due to molecular impacts is: 

dT/dt = (aE-Hj/mCp (1) 

The accommodation coefficient, a, is the fraction of the beam energy transferred 
to the crystal. The energy, E, carried by the beam, and the heat, H, dissipated 
by the crystal are in watts. The mass of the crystal is m in g, and Cp is its 
heat capacity in J/g°C. 

A drawing of the crystal is shown in Fig. 1. One side of the crystal is 
plated with the test material on which energy transfer measurements are to be 
made. If the plating is a conductor, it can also be used as one of the elec- 
trodes needed to drive the crystal in oscillation, otherwise, a conducting 
electrode must be first plated onto the crystal. The total amount of material 
that can be plated on a crystal before it significantly damps its oscillation 
is approximately 10,000 A for a light material such as aluminum and 1,000 A 
for a heavy material such as gold. Nearly the entire surface of the crystal is 
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Fig. 1 Quartz crystal temperature sensor 
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plated to permit measurements with beams having a circular cross section of up 
to 0.1 cm in diameter. 

Another electrode must also be plated on the backside of the crystal, to 
complete the driving circuit. It has a radius of 2 mm as measured from the 
center of the crystal so as to keep the plating capacity low and assure a good 
oscillating amplitude. 

The calibration tab is used to drive an electrical current through the 
plating so that the frequency change of the crystal can be determined when a 
known amount of power is dissipated at its surface. 


H = Ckiai/si + k2a2/S2)CT2 - T^) + (kjas/ssKTi - Tq) + a(T2^ - Tq4) (2) 

The heat conductivity of quartz, the plating, and wires holding the crystal are 
k^, k2, and k3 respectively in J/*C cm sec. The cross sectional area normal to 
the heat flow is ai , a2, and 33 in cm2. The conducting path length along the 
direction of heat flow is si, S2» and S3 in cm. The temperature at the center 
of the crystal where the beam impacts is T2. The temperature at the edge of the 
crystal is T^, and Tq is the ambient temperature. The Stefan-Boltzmann constant 
is o . 

For short bombardment periods 0, and the energy transferred to the 
crystal can be found from measuring its tenperature change. It was calculated 
that for a bombardment period T < 10 sec, H < ,03 aE and, 

aE = mcp dT/dt (3) 


to within 3 percent. 

Knowing the energy carried by the molecular beam, the accommodation co- 
efficient can also be found and: 

a = (mcp/E) dT/dt (4) 

Experimentally it was shown that the energy dissipated by the crystal over 
short periods is much less than the energy transferred from the beam for tempera- 
changes of less than one °C. By limiting the input of electrical power to 1 mW/cm2, 
the first 10 seconds of data is within a few percent of the calculated linear 
frequency increase . 

3.0 THEPMAL ANALYSIS OF THE ENERGY TRANSFER PROBE 

A thermal analysis was made on the energy transfer probe to optimize its 
operation. The probe utilizes an oscillating quartz crystal as the temperature 
sensor. Fig. 2 shows the probe. Four of these probes were flown on OGO-6. The 
probe was pointed into the upper atmosphere stream and when the shutter is opened, 
impacts of molecules transfer energy to the sensor crystal, raise its temperature, 
and cause its oscillating frequency to change. The change in frequency is pro- 
portional to the energy transferred. From the frequency change the accommodation 
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Fig. 2 OGO-6 Energy Transfer Probe 




coefficient of the surface plated on the crystal can be determined. The analysis 
includes the use of a BaF2 shutter so that the probe can be used in later work 
now being planned with neutral beams. 

Thermal analysis of the probe was accomplished by dividing it into a net- 
work of thermal nodes. The nodal network was solved by the solution of a set of 
linear first order differential equations. 

Design goals were: 

a. Optimize plating of the crystal sensor to reduce radiation losses. 

b. How best to coat instrument walls to reduce space temperature 
effects on the crystal. 

c. Effects of electronics leads on heat conduction. 

3.1 THEmL .ANALYSIS COMPUTER PROGRAM 

The analysis involved the use of several digital computer programs. The 
term "thermal model" refers to the basic physical parameters of a system which 
influence its tenperature within a given environment. Some of these parameters are: 

a. Mass 

b. Specific heat 

c. Thermal conductivity 

d. Thermal radiative surface properties (reflectance, transmittance) 

e . Geometries shape 

The probe was divided into parts called nodes. This is usually done on 
the basis of isothermal! ty. That is, those parts of the object which are con- 
sidered by the analyst to have a common temperature, due to high internal conduction 
paths, are lumped together. The degree of accuracy of the calculation procedures 
is very dependent on this process. In general, the more nodes considered, the 
more accurate the results . 

Temperature is a measure of the stored energy of an object. The process of 
computing tenperature as a function of time is an energy balance between the 
foiling factors : 

a. Heat energy lost or gained by conduction to or from another node 

b. Heat energy lost or gained by radiation to or from another node 

c. Heat energy gained by radiation from the radiative environment 

d. Heat gained through molecular energy accommodation to the surface 

e. Heat gained by the dissipation of electrical energy as heat within 

the node 

f. Heat lost by radiation to the environment 

When the above factors are summed, they produce a net gain or loss of energy 
to the node. Energy gained is stored in the mass of the node and results in 
positive change of tenperature. This change of tenperature is directly propor- 
tional to the product of the mass and specific heat of the node. Energy lost to 
the node results in a negative change of temperature with the same constant of 
' proportionality. 

The integration program is a logical procedure for evaluating the heat 
energy effect on a multinodal object. It can solve systems involving as many as 
200 nodes, considering 25 conduction paths and 25 radiation paths to each node. A 
maximum number of 400 conductive paths and 100 radiative paths is permitted in all. 
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The program is written for the greater part in the FORTRAN language, which 
permits easy interpretation or modification for special uses. 

3.2 NUMERICAL INTEGRATION 


The general equation for the storage of heat in mass is: 


This first order differential equation is replaced with the difference equa- 
tion, and with some rearrangement yields: 


Q^0 


( 6 ) 


The Eulerian technique of numerical integration yields the following equa- 
tion for the advancement of temperature at node (i) and time (0): 


Qt ® 

T(i, e * M - Tj,,9 . 


( 7 ) 


Temperature values on the interval 0 min to 0 max converge to the solution 
of equation (5) for small enough values of A0 . 

This numerical procedure is explicit, and convergence must be 
insured by check runs at smaller intervals of integration. However, if the 
results plot as smooth curves (stable solution) , the solution obtained is 

usually convergent. . 

The procedure of equation (7) is applied at every node based on 
heat flux conditions prevailing at the beginning of the interval. Tempera- 
tures for all the nodes and time are advanced to the end of the interval 
where heat flux calculation can be made based on the new set of temperatures 
based on the same recursion foimula, equation (7). 


3.3 HEAT FLUX CALCULATIONS 

The following heats are added algebraically to produce Qt (total 
heat) at node (i) . 

The dissipated heat, Qdig in the program, is an input value for 
each node considered. For those nodes in which no electrical energy is dis- 
sipated, the Qdis column is left blank. The value of Q^jis subject to an 
o£f-on switch in the environment table. System duty cycles may be simulated 
in this manner to turn the entire electrical system on or off at predetermined 

times . Qhis is entered in watts . ... 

Two physical surface properties of materials are important in com- 
puting the heat absorbed from the various energy sources mentioned above. 

These are the radiant energy absorptivity (a) and the thermal emissivity ( 0 ) 
of the node or its surface coating. The absorptivity measures the fraction 
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t 


o£ energy absorbed by the node (i) , as compared to that absorbed by a black 
body where the energy source has the spectrum of the sun. The thermal emis- 
sivity measures the fraction of energy emitted at node (i) as conrpared to 
that emitted by a black body of the same temperature. Where an energy source 
is in the same teirperature range as node (i) , the thermal emissivity (e) also 
measures the fraction of the source energy absorbed as conpared to that which 
would be absorbed by a black body. 

The tabular data is interpolated for the proper compute time by 
the following linear interpolation formulas: 


S 9 = ^ ~ 

^(i) ■ A9^ pi) 




- '’t 

Cl) 


^^Ci) 


0 0 -9t 

A0^ 


^Ci) 


0 0 -6t 

A0, 


0. + A0^ a 

Al(i) ■ ^^(i) ^ 


9^ + AO^. e^. 

^(i) -"^(i) 


Total environmental heat at time t is then ccmputed by: 


'(i) 


Ci)f^(i) ®(i)'^*^Ci) 


C8) 

C9) 

( 10 ) 


( 11 ) 


The heat which is exchanged by conduction is computed using: 



- hj)Sil 

hi) 



( 12 ) 


where m values of conductivity (K) conductive path cross section (A^) path 
length (1) , and the integers (a) and (b) are entered as input by the analyst 
for conductive paths numbered j =1 - m. The integers (a) and (b) denote 
the numbers of nodes (i) linked by the conductive path at (j) Block 2 which 
is called the ”Q Table" lists after the designator QC those conduction paths 
by number (j) which affect node (i) . The "Q Table" is developed by the 
program from conduction and radiation path input data. 

Heat which is exchanged by radiation between nodes is computed 

using : 


(k) 


= ^ ^ (k)\ 


(k) 



(13) 


where m values of script F ((? ) and the integers c and d are entered for 
each radiative pat k = 1 -* n in the same manner as for conductive paths. 
The designator QR proceeds the "Q Table" listing of radiative path numbers 
(k) for nodes (i) . The values of (:^) are obtained from the script (:?•) 
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program when an internal radios ity network is involved, 
to space may be handled by setting: 

^ (i] " ^(i) ®(i) 


External radiation 
(14) 


where F is the geometric view factor of node 
as a dtcnmy node of T = -46CPF. 

Qt is computed using: 

‘5'd) ■ 


(i) for space. Space is entered 


i i 

+ + ^QC (15) 

k' k 


That is, Qt is the value of electrical power dissipated at node (i) plus 
the environmental heat at (i) , plus the sum of the radiative heats at (i) , 
plus the sum of the conductive heats at (i) . 


3.4 QUASI STEADY STATE OPTION 

For systems which have large variations in the thermal mass 

of the various nodes, this option may prove of value. Nodes of 
extfmeiy small thermal mass tend to respond quickly to their thermal en- 
vironment. Where these light nodes are mixed with other nodes of greater 
mass, the integration step size (A0) for convergence is determined by the 
lighter nodes. This effect sometimes causes excessively long running time 
for the computer. To avoid this occurance, the analyst may sometimes wish 
to indicate some nodes as Quasi Steady State Nodes. The program then treats 
these nodes as though they had no mass at all so that they respond instan- 
taneously to their environment of the moment at each time step. The mathe- 
matical procedure used for accomplishing this effect is called Relaxation. 

The set of differential equations for the Q.S.S. nodes becomes a set of 
algebraic equations. The temperatures at each node are adjusted by a pro- 
cedure called "Regula Falsi" to drive the Q^^^ for each Q.S.S. node to zero. 

3.5 RADIANT EXCHANGE FACTOR PROGRAM 

Radiant losses from the Energy Transfer Probe are a function of 
the emittance and reflectance of the various surfaces which "see" the crystal 
(Fig. 3). In order to properly analyze the radiation losses to space, the 
effects of multiple reflection within the probe chamber must be considered. 
This analysis was performed using a radiant exchange factor program which is 
described below. 


3.6 


DEFINITION OF THE RADIATION FACTOR 


Radiant energy exchange between elements of a system of non-black 
bodies may be computed by the following equation: 




( 16 ) 


where , 
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Fig. 3 Theimal nodal diagram of energy transfer 


a = the Stefan-Boltzman constant 

Fjj = the ordinary view factor for black bodies 

Ai - the area of the i'th element 

Tj = the temperature of the i’th element 

Tj = the temperature of the j ’th element 

Fgij is a factor which accounts for the energy transferred between 
i and j where*^heat transfer is a function of multiple reflections between i 
and j themselves and other elements of the system. Fgj^j is a function of 
both the geometry of the system and the emissivity of tne various elements . 
The product of the geometric view factor Fjj and Fg^j is here called ij 
(called script F) . Equation (16) becomes 


3.7 



ij 


A, (T 



ASSUMPTIONS AND EQUATIONS 


(17) 


. The following procedure is adapted, with slight modification, from 
Gebhardt” . 

The total emission from a non-black surface may be represented by: 


qi . CjoT. Aj 


(18) 


where = the total hemispherical emissivity of object i. 

If j is then defined as the fraction of this total energy emitted 
by element i which is eventually absorbed by j , then the net interchange be- 
tween i and j can be expressed as: 


K . = B. e . - cT A.T. - 

^ij 1 ji 11 






(19) 


Since, for any geometric system, the following reciprocity relation must 
hold. 


B. . e -A. = B. . e -A. 

11 1 1 11 1 3 

By substitution of Equation (20) in (19), Equation (19) becomes: 


'‘ij = 'ij ''(^1 - T^) 


( 20 ) 


( 21 ) 


By comparison with Equation (17), it is apparent that: 



may be detemined in the following way. Consider The total 

fraction of energy leaving 1 and absorbed by j will be: 


®ij ' hj 'j * hi h hj * h2 ^ 


F, p_ B . 
In n nj 


( 22 ) 


The first term on the right side of Equation (22) represents energy trans- 
ferred directly, where Fij is the view factor of 1 for j and the total ab- 
sorptivity of j equals the total emissivity of j . 


e. =a. by Kirchoff’s law 

Here all surfaces are considered "grey bodies" idiich have the same value of 
spectral absorptivity at all wavelengths and which also are isotropic emitters. 

The third term on the right hand side of Equation (22) is typical 
of all other terms where F ^2 is the view factor of 1 for 2, P2 is the reflec- 
tivity of 2, and B 2 j is as previously defined. Note that the assumption is 
made that element 2 reflects energy in the same manner as it emits energy 
(isotropic emission). That is, that there is no preferential direction for 
reflection (diffuse reflection) . Since an equation similar to (22) may be 
written for each element of the system, a linear system of equation is de- 
fined with some rearrangement of equation (22) . 


(Fll Pj - 1)B.^ ^ 12^2 * •••* ^In ‘’n " ' ^Ij 

^21 ^l®lj ^^22 ^2 ‘ ^^®2j ^2n ^n \j “ ' ^2j 


F T p B . • + 
nl ^n i;i 


* ^^nnP n ^^®nj ' ^nj ® j 


The coefficients of the unloiowns may be expressed in matrix fashion by: 

= Fi^ Pj - 6 (23) 


where 6 is the Kronecker delta. 

The n values of Bij for each j are a function of the right hand 
side column vector ^ij ■ 'fij €j. Since, however, n^ values of B are re- 
quired for each j = l"^-* n for which the M matrix does not change, then in 
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matrix fom. 


B = M"^a 


(24) 


where B and a are square matrices of order n. 


Pi = 



(25) 


for opaque grey surfaces. 

In order to obtain a square array of j values for use in Equation 
(17), each row of the B matrix may be multiplied by the emissivity corres- 
ponding to the elements i = 1 -* n. 

3.8 . DETAILED THERMAL MODEL 

Thermal models were created for both 10 MHz and 15 MHz crystals. 

A general schematic of the nodal system for the crystal is given in Fig. 4. 
This schematic applies to the test instrument using 15 Wz crystals. The 
model using 10 MHz crystals, differs somewhat from this. The length of the 
steel crystal supports were longer for the 10 MHz crystal. The third lead 
was not well defined for the 10 MHz crystal. The surfaces of the crystals 
were divided into 21 nodes, conductive paths from the crystal edges to the 
mounting were divided as shown. 

Fig, 5 is a plot of the comparison between laboratory tests on 
the 10 MHz crystal and a computer simulation made of the same situation. 
Power dissipation was 760 iJi watts. The center curve of this figure shows 
the laboratory measurements of frequency response when electrical energy 
was dissipated in the gold plating on the crystal face. A third lead was 
used to conduct power to the crystal , Unfortunately the heat leak along 
this wire is not well defined since an alligator clip connector was used. 

For this reason the model was run assuming 

a. no heat leak along the third lead (top curve), and 

b. assuming a leak to match that of the support leads. 

Fran this figure it can be seen that the initial slopes of all three curves 
are well matched. The slight delay in the laboratory curve is due to a slow 
collection rate for the data acquisition system. It is readily apparent 
that adjustment of the third lead leak can be used to align the curves al- 
most perfectly. 

Fig. 6 compares time temperature response curves for both the 
37 node model and the simplified 12 node model of the same device. 

Fig. 7 shows the third cycle of computer predicted instrument 
response to a chopped molecular beam, where 300 ijW/cm^ are absorbed by the 
crystal surface. In simulation the instrument was run to steady state with 
the window open. Then the beam was imposed and the chopper window was al- 
lowed to follow its normal cycle. It can be seen that in three cycles the 
instrument approaches a cyclic equilibrium. Computer runs used during the 
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Fig. 4 Thermal nodal diagram of sensor crystal 
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study were 

a. Run to shutter open equilibrium, no beam 

b. First cycle, shutter open, beam on 

c. First cycle, shutter closed 

d. Second cycle, shutter open 

e. Second cycle, shutter closed 

£. Third cycle, shutter open 

g. Third cycle, shutter closed 

4.0 EXPERIMENTAL MEASUREMENTS 

With the thermal analysis o£ the probe completed, experiments were 
beg\m on the bombardment of metals in beams of noble gases. Simultaneously, 
measurements o£ sputtering yields and kinetic energy transfer rates were 
made at energies between 50 eV and 600 eV. This energy range is essentially 
unexplored because most measurements are being done below 10 eV for applica- 
tion to hypersonic flow problems. 

Impacts studies between 50 eV and 600 eV are important because 
this is a transition region where sputtering yields and energy transfer rates 
increase rapidly with the particle impact energy. This follows because at 
the low end of the range interactions take place with the first monolayer of 
target atoms and bombarding particles tend to be reflected with little energy 
transfer. At the high end of the range, the bombarding particles penetrated 
a few monolayers into the target and become momentarily trapped and energy 
transfer is nearly complete. 

Measurements were made at normal incidence on polycrystalline A1 
and Au in Ne'*’, and Ar"*" beams. The surfaces were plated on a 1/2 inch dia- 
meter mass and temperature sensitive optically-polished quartz crystals. 

Mass changes of 5 x 10’^^ g/cm^ and tenperature changes of 6 x 10“^ C were 
detected. 

The results for A1 in Ne'*’ and Ar"'' are shown in Figs. 8 and 9. 
During the bombardment period of 200 seconds , target temperatures increased 
by as much as 0.7 ®C above the ambient. 

Results for the bombardment of Au in Ne'*’ and Ar"*” are shown in 
Figs. 10 and 11. 

Data from Figs. 8-11 were used to obtain the energy accommodation 
and sputtering yields for the ion-target combinations and are shown in Figs. 
12 and 13. 


4.1 ENERGY ACCOMODATION 

For energies below 600 eV, a classical approach to the computation 
of the accommodation coefficient to predict experimental results can be used 
because the de Brolie wavelength is much smaller than the lattice spacing. 
For the data given here, a Bohr^ potential with hard-core approximation can 
be used to describe the collision that occurs between the impacting particle 
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Fig. 8 Temperature rise of A1 in a 10 ^A/cm^ Ne^ beam at different bombardment energies 
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Fig. 11 Temperature rise of Au in a 10 uA/cm^ Ar'*’ beam at different bombardment 
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and the target atom. 

The transfer of energy is bounded within limits of from zero to 
4 Emxm 2 /(m]^ + m 2 )^ where the impact energy is E, the incident particle mass 
is and the target atom mass is m 2 . Taking zero as the lower limit, the 
resulting average energy transfer is 2 m^m^E/Cmi + ^ 2 )^. By dividing the 
expression by E, the accommodation coefficient, a, for the surface is ob- 
tained. Such an analysis assumes only one collision with the target solid 
and subsequent return of the incident particle to the gas. Although near- 
ly all of the work reported in the literature has been at energies below 
10 eV, a short discussion of this work is valuable here to see how accommo- 
dation measurements at low energies compare with our initial results at 
higher energies. 

The simplest theory of how particles transfer energy to a surface 
was published by Baule® in 1914. He applied classical theory and assumed 
a gas model with a Maxwellian distribution each of which make a single col- 
lision with the target solid before returning to the gas, and determined 
the energy accommodation coefficient, a, to be: 


a 


2 

(mj^ + m2) ^ 


(24) 


Experimental work with low temperature gases has shown marked dis- 
agreement with Baule's approximations. Consequently, workers in the field 
have since become highly sophisticated to bring experiment and theory into 
closer agreement by considering crystal responses to atomic impact and the 
attractive and repulsive portions of scattering potentials in addition to 
the simple hard sphere repulsion. q 

Theoretical work done by Goodman^ on the response of a three-di- 
mensional crystal to a noimally incident atom fits our experimental work to 
SOTe degree. He finds the high energy limit for the accommodation coeffi- 
cient as 4 mim 2 / (m^ + m 2 )‘^. This result is not too surprising because only 
head-on collisions were considered with the collision time much less than 
the natural period of the crystal lattice. The last condition is tantamount 
to neglecting the lattice coupling insofar as energy transfer during colli- 
sion is concerned, 

Trilling^ also considered collisions along the line of centers 
but allowed all angles of incidence. He computed the accommodation coeffi- 
cient at high energies a « 1.78 mi/m 2 provided mj^/m 2 » 1. This result can 
be understood qualitatively from GooOTan's result by averaging over the solid 
angles of incident direction. 

Values in Figs. 12 and 13 reach the maximum energy accommodation 
given by Goodman at 600 eV. We believe, though, the large transfer of energy 
resulted from particle trapping and not from neglecting of lattice coupling. 
Investigations are now being made on the penetration range of impacting par- 
ticles as a function of impact energy to determine at what energy the impact- 
ing particles are effectively trapped. 
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1 .0 INTRODUCTION 

At apogee OGO-6 periodically passed through belts of high- 
energy particles. The particles are composed mainly of electrons 
and protons. OGO-6 surfaces will be damaged by these particle impacts 
that will penetrate the surface and displace atoms leading to the dis- 
coloration of optical coatings . 

When a surface is bombarded by high-energy particles, it will 
first show a net mass gain from the trapping of particles and then a net 
mass loss from erosion as successive layers are sputtered by the high- 
energy impacts. The mass rate of change of the surface, dM/dt, can 
be expressed as 


dM/dt = Of m^i- p (t)m^ - 

The absorption coefficient a is the fraction of ions that pene- 
trate and become trapped in the target. The re— emission coefficient 
^ (t) is the number of trapped particles re-emitted as succeeding 
layers of target atoms are sputtered. The flux is in particles/cm^-sec . 
The number of surface atoms eroded by particle impact is p . The 
mass of the bombarding particle and of a surface atom are, respec- 
tively, m-j and mg. At t = 0, p (t) = 0, and as t becomes large, 
p(t) _ Of i . Once P(t) = of i the target will lose mass at a constant rate 
due to sputtering. 

To penetrate beyond the first layer of surface atoms, the 
inequality 


CT n 


2/3 


^ 1 


( 2 ) 


will hold, otherwise the particle will probably be reflected. With n 
the atomic number density of the surface, n^/® represents the area 
number density seen by the incident particle with a total binary cross 
section, cr . This inequality is satisfied for protons and alpha particles. 

Although the cross sections at these low energies are not known 
accurately, Bohr^ has shown a reasonable cross section to be 
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which is replaced by 
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n R 


( 4 ) 
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whenever cr >> /(Z^ + Z 2 ). At the lowest energies the binary 

collision cross sect?on approaches the gas kinetic value. R is the distance 
of closest approach obtained as the solution of 




= 2Z Z^E^a expr-RCZ^"^^ + Z^^®)^^^/a ] (5) 

2^ 1 2 R o ^ 1 2 o-* ^ 


Subscripts 1 and 2 designate incident and surface particles, respec- 
tively, with E the kinetic energy, m the mass, A the atomic number, 

Sq the Bohr radius and Ep^ the Rydberg energy. 

To obtain the distribution of the penetrating ions it is necessary 
to consider their motion in the surface. Although the ions are probably 
neutralized upon entering the surface this is expected to be of no con- 
sequency since hard sphere type collisions are assumed to occur. For 
protons and alpha particles, m 2 /m^ v, 1 , large angle scattering occurs 
at each collision giving an erratic rectilinear path which resembles a 
diffusion process. In fact, the slowing down process is treated like that 
of neutrons where the particle slowing down density for a plane flux, i, is 


-X 


q = 


4R 


(4TT R) 


1/2 


( 6 ) 


p 

where R is the "Fermi age" . The density, q, represents the number 
of particles that slow down past an energy, E, at a given depth, x, in 
the surface* If rr\ 2 /^^ > ^ f the angular deflection of the incident par- 
ticle from its direction of travel is small and the particle straggles 
about the penetration depth with an approximate Gaussian distribution. 
Because the notion of range is complex, it is customary to 
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define a linear, a vector and a projected range. The linear range 
is the total rectilinear path traversed by the particle to the point of 
stopping while the vector range is the vector distance from point of 
entry into the surface to the end point. Projection of the vector dis- 
tance along a line normal to the surface gives the projected range. 
An approximate relation between the average linear range R|_ and 
the average projected range Rp (or mean range) was developed by 
Linhard and Scharff® as 
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1 "^2 

Rd (1 

P ^ 3 m 
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( 7 ) 


for an inverse square potential. For mgm^ < 1 » approximating (7) is 
quite accurate but, for m 2 /m^ ^ 1 the approximation has application to 
a portion of the low energy region only.'^ 

Scattering near the end point of the particle travel is probably 
a man-body collision problem and not the binary one assumed to obtain 
the distributions discussed in the previous section. If relation (2) is 
used as a necessary criterion for halting the motion of the incident 
particle, the cross section must be approximately the gas kinetic value, 

^ 5 X 10“^® crn2; cross sections of these magnitudes correspong to 

energies below 1 .0 eV for protons and alpha particles. At these low 
energies the de BroUe wave length is comparable to the lattice separation 
indicating that more than one target atom contributes to a collision. 

Andersen and Sigmund® have presented a model for non-binary 
scattering which is of importance at much higher energies than indicated 
by (2). It recognized the fact that, with such small mean free paths, 
correlations between successive collisions cannot be neglected in the 
low energy region. The model assumed involves a set of crystal lattice 
rows which are resolved into stacked symmetric rings of atoms. The 
incident particle moves along the symmetry axis of these rings with the 
projectile-ring interaction as the basic event with the ring presenting 
an effective potential field to the incident atom. 

If the incident particle is aligned along the symmetry axis, close 
collisions become improbable because the field keeps the particle near 
the symmetry axis to produce channeling. For the interpretation of 
penetration distributions known in the literature as ’’supertails" it is 
important to investigate the low energy collisions to determine the 
density of vacancies, interstitials and self -interstitials created. From 
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this, estimates of further projectile penetration by diffusion after 
the slowing down process can be made. Using a Born-Mayer potential 
for each lattice atom, a new reflecting (or stopping) energy replaces 
the one determined through (2) and is found to be^ 

R „ p Ae“^^^ (8) 

refl 

whenever m^ » m 2 * Constants p and L define the number of particles 
in the ring and radius of the ring element} A and a are the potential 
parameters . Reflection energies given by (8) are somewhat higher 
than those found through (2); for alpha particles in heavy materials, 
typical estimates are 5-10 eV. 


2.0 DIFFUSION OF TRAPPED PARTICLES 

Although channeling gives a plausible explanation of penetration 
of about one micron for these low energy particles, the "supertails" 
observed® extend to depths of several microns. A possible explanation 
to these deep penetrations is the steady state diffusion approach of 
Sparks.^ He assumes the diffusion of vacancies and self-interstitials 
created by the Incident particle is unaffected by the relatively small 
number of implanted ions. With the bombarding particles forming 
substitutionals through the annihilation of vacancies two general forms 
of the integral penetration distributions, hgCx), result which may be 
expressed as 


a exp[-k(x-i )]} X > (9) 


and 


.2+r 


n_(x) 


(x -l-V d) 




( 10 ) 


The constant, d, depends upon the vacancy density interstitial density 
and the diffusion coefficient, i , is located approximately at the edge of 
the stopping region. Distribution (10) normally occurs if vacancies 
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and self-interstitials are created and annihilated in pairs. For 
distribution (9) to obtain the vacancies, the self-interstitials must 
be created at different rates which means that most of the inter- 
stitials must be of the bombarding particle type . 


3 . 0 BOMBARDMEN T OF A LUMIN UM 

Laboratory measurements were made on the bombardment 
of aluminum surfaces by protons to determine the saturation concen- 
tration of protons in this material commonly used on OGO-6 . The 
function, q, given by (6) was used to calculate the distribution . The 
distribution is for the case m-j > mg. It has a maximum near the 
surface because of large angle scattering that occurs for the case 

in which the incident particle mass m^ is less than the target atom 
10 ' 

mass mg . 

■^e mean penetration depth Rp of the particles needed to 
scale the distribution function can be found from the p -e relation- 
ship of Lindhard and Scharff , ^ ^ which relates the range-energy of 
particles in solids: 


P 


e 


= R 


166 


m (1 + 3Ap ^^2/3 ^2/3^ ^ 


33 




( 11 ) 


( 12 ) 


where the subcripts 1 and 2 refer to the incident ion and the target 
atom, and is the energy of impact in keV. 

Davie's group has shown that the P - e relationship is accurate 
to about 20% for € >10. Values for p , knowing s , have been 
plotted by them for 0< e <10 and 0< p <100.^^ 

The absorption coefficient a was found, experimentally, to 
be about 0.01. 

The saturation concentration of trapped hydrogen in aluminum 
is shown in Fig. 1 . 
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Fig. 1 Saturation concentration of hydrogen 
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ABSTRACT 

A quartz crystal microbalance is of limited value in monitoring 
surface contamination on satellites or in space simulation cham- 
bers because it operates several degrees above ambient tem- 
peratures. The amount of contamination adsorbed on a surface 
is highly temperature dependent and the higher temperature of 
the microbalance will significantly reduce the amount of con- 
tamination it adsorbs. Generally, a quartz crystal microbalance 
will indicate a lower level of contamination than the amount that 
is actually present. To overcome this problem, a thermoelec— 
trically-cooled quartz crystal microbalance has been developed 
to monitor surface contamination as a function of temperature. 

INTRODUCTION 

The use of a quartz crystal microbalance (QCM) is a well 
established method for weighing thin films of solid materials, down to 
a fraction of an Angstrom, that are deposited"* or removed^ from a 
surface. Solid films couple strongly into the oscillating QCM and its 
frequency change is proportional to the mass loading. Direct applica- 
tion of a QCM to monitor a wide range of contamination in space simu- 
lation chambers has proven to be a more difficult task because a QCM 
operates several degrees above ambient temperatures . Gaseous con- 
tamination in equilibrium with solid materials in the chamber is not 
readily adsorbed on a QCM because of its higher temperatures . Con- 
tradictory results are often obtained where a passive optical system 
shows contamination while a QCM operating nearby indicates little is 
present. 

As the study of surface contamination becomes more fundamen- 
tal in nature, a QCM specifically designed to monitor contamination is 

* This work supported under NASA Contracts NASS— 27879 and 

NAS5-1 1163. For presentation at the AIAA— NASA/ASTM/IES 7th 
Space Simulation Conference, November 12-14, 1973, Los Angeles, 
California 
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needed. Not only must the heat generated by the QCM be removed, 
but a convenient method found to automatically control its operating 
temperature, so that, surface contamination studies can be made 
over temperature extremes common to spacecraft. To permit in-situ 
reflectivity measurements, the crystals should be optically polished. 

With these objectives in mind, a thermoelectrically-cooled quartz 
crystal microbalance has been developed for monitoring contamination 
on optical surfaces as a function of temperature. 

QCM LIMITATIONS 

A QCM is an active quartz crystal oscillator and power dissipation 
raises its temperature several degrees above the ambient. Because 
the amount of contamination adsorbed on a surface is highly temperature 
dependent, a passive surface at ambient temperature will adsorb a 
number of contamination monolayers until it reaches an equilibrium 
with sources in the chamber while a QCM will adsorb a smaller number 
because of the different equilibrium condition. In some cases, a heavily 
driven QCM will operate about 10“C above ambient. At this elevated 
temperature no measureable amount of contamination will be observed. 

The amount of contamination adsorbed on a surface^is dependent 
on the residence time, T , of the contamination molecules . 

T = T ^ exp ( A E/RT) 

where t o vibrational period of the contamination lattice, A E 

the desorption activation energy, R the gas constant and T the absolute 
temperature. Griffith‘S has shown that the desorption rate of contamina- 
tion with desorption activation energies of less than 25,000 cal/g— mol 
(oils and epoxies) is highly dependent on temperature. A 10°C tempera- 
ture rise can result in nearly a ten fold increase in the desorption rate. 
Use of a QCM to monitor contamination will consistently give lower 
readings than the actual level present because of its higher temperature. 

The strong effect of temperature on the adsorption and desorption 
of contamination is shown in Fig. 1 . Here long-term OGO-6 measure- 
ments are shown correlating QCM contamination loading to the eclipse 
period of the satellite. As has been previously reported, the primary 
source of contamination on the satellite was the solar panels baking out 

5 

in the sun . 

OGO-6 was inserted into a polar orbit with its orbit plane normal 
to the earth-sun line. This orbit was chosen so that the satellite would 
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be in full sunlight for the maximum period after launch. During the 
four week period after launch, the amount of contamination adsorbed 
onto the QCM steadily increased until the first week in July when it 
abruptly decreased and finally in August the QCM showed a net loss 
of contamination and became a contamination source. By correlating 
the eclipse periods of the satellite to the QCM measurements, the 
reason for the fluctuation in the contamination adsorption and desorp- 
tion rates became apparent. 

In full sunlight the temperature of the solar panels was 72“ C 
and a wide range of high and low volatile contamination outgassed from 
the panels onto the QCM. During maximum eclipse, when OGO-6 was 
in the earth shadow 30% of the time, the average temperature of the 
solar panels dropped to 60° C and the contamination flux from the panels 
decreased significantly. The QCM lost contamination during the eclipse 
period because tine lower flux rate from the panels did not balance out 
contamination desorbing from the QCM. The QCM desorbed contamina- 
tion into space and onto adjacent surfaces that were at a lower tempera- 
ture than the QCM . 

Although the QCM indicated a loss of contamination during eclipse, 
contamination was present. Reber's Neutral Mass Spectrometer, located 
next to the QCM, showed that a strong flux of contamination continued to 
outgas from the panels. 

THERMOELECTRICALLY-COOLED QCM 

In order to conduct contamination experiments as a function of 
temperature, a new instrument called a thermoelectrically-cooled QCM 
(TQCM) was developed. A thermoelectric device was picked to control 
temperature because it offers several advantages over other methods. 

It uses only a series of solid-state bismuth telluride junctions, through 
which electrical current is passed, to pump heat to or from a load by 
the Peltier effect. It is small in size and can be remotely operated 
with a signal pair of electrical leads. It has no moving parts and is 
highly reliable. There are no requirements for pumping a refrigerant 
or for supplying a coolant, such as, liquid nitrogen. These features 
result in greatly reduced operating and maintenance costs in controlling 
the temperature of a QCM with a thermoelectric device. 

The design of the TQCM has been previously reported® and only 
its operation will be described here. The TQCM instrumentation opera- 
ting under ambient conditions is shown in Fig. 2. The TQCM Controller 
operates on 115 Vac, 60 Hz at 0.25 A and provides the various voltage 
outputs to operate the crystal oscillators, temperature bridge. 
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Fig. 2 TQCM Operating at +15°C under Ambient Room Conditions 


thermoelectric device and readouts of temperature and frequency. 

The TQCM Electronics, Heat Sink, and Sensor are shown in 
Fig. 3. The TQCM uses a two-stage thermoelectric device to auto- 
matically control the sensor crystal temperature between -50° C and 
+ 100°C to i 1 °C in vacuum. For operation over this range, the heat 
sink temperature is to be maintained below +40°C. No problem will 
be encountered in maintaining the heat sink below +40°C by mounting 
it on a 0.6 cm thick metal bracket. Heat generated by the TQCM is 
readily dissipated by the bracket at ambient room temperatures . 

The TQCM has been designed with an extended temperature 
bridge circuit and reserve power and will reach a lower temperature 
limit of -59°C if its heat sink is maintained below +25°C, To maintain 
a +25°C temperature, the heat sink is to be mounted on a 1.2 cm thick 
metal bracket capable of removing a maximum 2 . 8 W when the TQCM 
is at -59“C. The TQCM cool-down time in vacuum is shown in Fig. 4. 

The TQCM operating temperature is set by positioning the 
thumb wheel switch on the Controller. A 3 1/2 digit panel meter is 
provided in the Controller for direct readout of temperature. 

A temperature output is provided for remote readout or recording . 
The TQCM temperature sensor is a precision platinum resistance ther- 
mometer linear to -.5%. At 100“C the voltage output is 1 ,00 Vdc. At 
-59°C the voltage output is - 0,59 Vdc. The TQCM operating temperature 
in °C equals 100 times the voltage output of the Controller. 

A matched pair of precision 10-MHz quartz crystals is used to 
measure mass loading. The crystals are designated as a sensor and 
reference crystal . The crystals are optically-polished and plated with 
Al. The sensor crystal is coated with magnesium fluoride for in-situ 
reflectivity measurements while contamination is collecting on its sur- 
face. The crystals can be changed by unloosening two set screws in the 
TQCM Sensor. 

The output frequency of the TQCM is the beat frequency between 
the two oscillating crystals. The beat frequency effectively eliminates 
frequency changes caused by ambient temperature variations. By care- 
fully matching sets of crystals in vacuum, frequency change with tem- 
perature of less than + 50 Hz between -59°C and +100°C is attained. 
Because only the sensor crystal sees the contamination flux, the TQCM 
output frequency will increase with mass loading . The crystal sensor is 
optically polished and full plated on one side . This technique produces 
a more active crystal whose mass sensitivity is greater by about 20% 
than semi— polished crystals. The TQCM mass sensitivity, m, is 

-9 2 

m = 3.5 X 10 g/cm -Hz 
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Fig. 3 TQCM Flectronics, Meat Sink and Sensor 




A frequency output is provided in the TQCM Controller to 
measure the frequency change of the crystal sensor produced by mass 
loading. Frequency is measured by a counter, such as, a Hewlett 
Packard Model HP5321B, 

The TQCM can be set at any particular temperature over its 
operating range. If the TQCM is set at ambient, the heat generated 
in the oscillating crystal will be removed so that contamination loading 
to passive optical systems can be monitored. By increasing tempera- 
ture to about +100^C the crystal can be cleaned. By periodically 
dropping the temperature of the TQCM in fixed steps, the amount of 
surface contamination for different equilibrium temperatures can be 
determined and the background level of contamination in a space 
simulation chamber monitored with time. 

The TQCM can be used to generate a calibrated contamination 
flux so that the contamination sticking coefficients of surfaces can be 
measured. A TQCM calibrated source is generated as follows. Con- 
tamination from an uncalibrated source is directed at a TQCM cooled 
to -SO'C. Contamination is allowed to freeze out on the TQCM. The 
uncalibrated source is then turned off and the vacuum system allowed 
to pump down. The TQCM now becomes the calibrated source by 
simply raising its temperature so that the contamination will desorb 
at the desired rate. By monitoring the TQCM frequency increase 
with time, the contamination mass flow rate can be accurately de- 
termined. 

The TQCM can also be used to calibrate the mass loading of 
viscous films on the quartz crystal sensor that do not couple well into 
the oscillating mass of the crystal. As a contamination film grows to 
several hundred monolayers, its top layers will slip relative to the 
oscillating crystal . The simple relationship between mass loading and 
frequency change will not hold true. Calibration is achieved by drop- 
ping the sensor temperature to freeze the contamination and provide 
rigid mass coupling to the crystal. 

CONTAMINATION MEASUREMENTS 

Studies of background contamination in a vacuum chamber were 
made to show the capability of the TQCM to measure surface contamina- 
tion under equilibrium conditions. A 1 50 l/sec Vac-Ion pump was used 
to evacuate the chamber to the 10”® Torr range. Just prior to the 
measurements, the chamber was contaminated with Welsh Duo-Seal 
roughing pump oil that raised its pressure into the 10"^ Torr range. 
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The net mass gain or loss for an aluminum surface from the 
adsorption and desorption of background contamination present in the 
chamber at 7 x 10“^ Torr is shown in Fig. 5. Below -30 °C adsorption 
predominated and contamination collected on the surface. Above 
-30°C desorption predominated. The break in the curve at about 
-10®C where the rate of desorption increased abruptly is significant. 

The break shows a phase change occuring in the contamination 
from a solid to a liquid. The phase change caused a large increase 
in its vapor pressure and the higher desorption rate. We were unable 
to make a measurement at +20°C because the desorption rate became 
so large that by the time the surface reached +20° C the contamination 
was nearly completely desorbed. 

Fig. 6 shows the desorption rate of roughing pump oil as a 
function of the number of monolayers on the surface for a contamina- 
tion background pressure of 4 x 10”^ Torr. When the number of 
contamination layers reaches 200, the increase in the desorption 
rate starts to level off indicating the bulk properties of the oil pre- 
dominate and the properties of the surface have little effect on the 
desorption. At 0°C the leveling off of the desorption rate is quite 
evident. 

The desorption rate decreases rapidly as the number of 
layers present is reduced because of the higher bonding energy holding 
contamination layers near to the surface. The higher bonding energy 
results from contamination molecules filling cracks and crevices in 
the polycrystal line aluminum surface. 

The number of monolayers of contamination on a surface in 
equilibrium with the contamination on the walls of the vacuum chamber 
as a function of temperature is shown in Fig. 7. The ambient tem- 
perature of the chamber was 23°C. The measurements were made 
by raising the TQCM temperature to +80°C to bake off surface con- 
tamination and then dropping its temperature to -50°C to adsorb a 
film of contamination. The TQCM temperature was then increased 
to a particular equilibrium temperature shown in the figure to detei — 
mine the number of contamination monolayers. From the figure it 
can be seen that there was always oil contamination on the surface 
even at 4 x 10“^ Torr for temperatures below +30°C. It would have 
been impossible to make the measurements with a QCM because there 
are no net mass changes at equilibrium . 

CONCLUSIONS 

Most contamination is adsorbed on a surface in the gaseous 
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Fig.5 Contamination adsorption and desorption rates for an Al surface 
with temperature in a 7xiO"^Torr vacuum 
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Desorption rate (lO'S/cm^sec) 



Fig. 6 Contamination desorption rate as a function of monolayers on the surface 
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Monolayers of Contomination 



Fig. 7 Number of contamination monolayers on a surface at equilibrium 
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or liquid state and the adsorption rate is highly temperature dependent. 
A QCM is of limited value in monitoring contamination because it 
operates several degrees above ambient temperatures and does not 
readily adsorb contamination. For accurate measurements, a tem- 
perature controlled QCM should be used. 

If a QCM is used to monitor contamination, measurements 
should be made from a directed source. The source temperature 
should be much higher than the QCM operating temperature to insure 
that the contamination sticking coefficient is greater than zero. 
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